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The Plateau-Rayleigh instability of a liquid column underlies a variety of fascinating phenomena 
that can be observed in everyday life. In contrast to the case of a free liquid cylinder, describing 
the evolution of a liquid layer on a solid fibre requires consideration of the solid-liquid interface. 
In this article, we revisit the Plateau-Rayleigh Instability of a liquid coating a fibre by varying 
the hydrodynamic boundary condition at the fibre-liquid interface, from no-slip to slip. While 
the wavelength is not sensitive to the solid-liquid interface, we find that the growth rate of the 
undulations strongly depends on the hydrodynamic boundary condition. The experiments are in 
excellent agreement with a new thin film theory incorporating slip, thus providing an original, 
quantitative and robust tool to measure slip lengths. 


Glistening pearls of water on a spider’s web [1], or the 
breakup of a cylindrical jet of water into droplets, are fa¬ 
miliar manifestations of the Plateau-Rayleigh instability 
(PRI) [2, 3]. By evolving into droplets, the surface area 
of the liquid, and consequently the surface energy, are re¬ 
duced. This instability also acts for a liquid film coating 
a solid fibre [4-12], though the flow boundary condition 
at the solid-liquid interface provides additional complex¬ 
ity to the system. While the breakup of a homogeneous 
him into droplets on a hbre may be a nuisance in coating 
technologies of e.g. wires and optical hbres, this funda¬ 
mental instability turns out to be very useful: take for 
example water collection through fog harvesting [13, 14], 
a biomimetic approach that is perfected in nature by the 
spider’s web [1]. Hence, taking advantage of the PRI on 
a hbre, enables gaining physical insight into the solid- 
liquid boundary condition. Indeed, understanding the 
break-down of the no-slip boundary condition [15] is of 
major interest to the scientihc and industrial communi¬ 
ties as it has practical implications in areas that involve 
small-scale huid systems, such as lab-on-a-chip devices, 
hows in porous media as well as biological hows to name 
a few. 

The classical case of the breakup of a laminar free 
howing liquid stream is dependent on material proper¬ 
ties and geometry [16]. While there are some studies of 
this instability in geometries other than free cylindrical 
liquid hows [17-23], the role of the boundary condition 
at the solid-liquid interface on the evolution of the PRI 
is less well understood [24, 25]. It has been shown that 
the appearance of a dewetting liquid rim undergoing a 
Plateau-Rayleigh-type instability is inhuenced by the hy¬ 
drodynamic boundary condition at the solid-liquid inter¬ 
face [26]; however, a general quantitative match between 
the growth dynamics of the PRI in experiments and an¬ 
alytical theory involving hydrodynamic slip has thus far 
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FIG. 1. Plateau-Rayleigh instability on a fibre. 

(a) Schematic and (b) optical micrographs illustrating the 
PRI for a liquid polystyrene film on a glass fibre. At t = 0, 
the polystyrene film on the no-slip fibre has a thickness 
eo = 13.2 ± 1 fim and the glass fibre radius is a = 9.6 ± 1 pm. 
The width of the optical images is 560 pm. 


not been achieved. 

Characterisation of slip boundary conditions at in¬ 
terfaces and their controlling parameters has been ac¬ 
tively investigated in the literature [27-31]. The classical 
boundary condition assumes no slip at the solid-liquid in¬ 
terface. That is, for simple shear flow in the x direction 
along a solid interface placed at z = 0, the tangential 
flow velocity, v„c(z), vanishes at the solid-liquid interface: 
^^a:(0) = 0. However, as already noted by Navier [32], 
there is no fundamental principle requiring Ua;(0) = 0, 
and one can also have hydrodynamic slippage at this in¬ 
terface, as defined by the slip length: b = [vx/dzV,c]\z=[), 
where 6 = 0 corresponds to the classical no-slip case. 

In this article, we explicitly address the effect of a vary¬ 
ing boundary condition, from no-slip to slip at the solid- 
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liquid interface, on the PRI of a viscous liquid layer on a 
solid fibre. We find that the growth rate of the instability 
is strongly affected by the solid-liquid interface, with a 
faster breakup into droplets for a slip boundary condi¬ 
tion compared to an equivalent sample with no-slip. In 
contrast, the wavelength A* of the fastest growing mode 
is not sensitive to the solid-liquid interface. The linear 
stability analysis of a newly developed thin film equation 
incorporating slip is in excellent agreement with our data. 
The theory is valid for all Newtonian liquids and enables 
the precise determination of the capillary velocity 7 / 77 , 
and most importantly of the slip length b. 


RESULTS 


The experimental approach. An entangled 
polystyrene (PS, 78 kg/mol) film with homogeneous 
thickness eg (5 to 93 /xm), is coated onto a fibre with ra¬ 
dius a (10 to 25 Tim), resulting in a PS coated fibre with 
radius ho = a + eg, as schematically shown in Fig. 1(a). 
Glass fibres provide a simple no-slip boundary condition 
[36]. In contrast, a slip interface results from coating the 
entangled PS film onto a glass fibre pre-coated with a 
nanometric thin amorphous fluoropolymer (AF2400, 14 
± 1 nm) [33]. The fluoropolymer coating on glass was 
used because it is well established that PS, above a crit¬ 
ical molecular weight (M^ ~ 35kg/mol), exhibits signifi¬ 
cant hydrodynamic slip at this solid-liquid interface [33]. 
Henceforth, we will refer to these as the “no-slip” and 
“slip” fibres. 

All samples were prepared and stored at room tem¬ 
perature, well below the PS glass transition temperature 
(Tg ~ 100 °C), thereby ensuring that there is no flow in 
the PS him prior to the start of the experiment. Prior 
to each experiment, the PS coated hbre was measured 
with optical microscopy as shown in the t = 0 image 
of Fig. 1(b). To initiate the experiments, samples were 
annealed in ambient atmosphere at 180 °C ~ well above 
Tg - which causes the PRI to develop. The evolution of 
the surface prohle was recorded with optical microscopy. 
Figure 1(b) shows a typical evolution for a PS him on a 
no-slip hbre. 

Above Tg, the liquid PS him becomes unstable, which 
causes variations in the local axisymmetric surface pro¬ 
hle, h{x,t) = hi) + C,{x,t), over the axial coordinate x 
and time t. The amplitude <^(x,t) of the undulations 
grows with time and hnally results in a droplet pattern 
displaying a uniform wavelength. A*. By measuring the 
spatial variation of C,{x, t) in the initial development, and 
locating the maxima, the PRI wavelength of each sample 
was determined (see Fig. 1(a)). Typically, four or hve 
wavelengths were averaged per sample. In addition, by 
measuring the temporal change in radius of an individ¬ 
ual bulge, we gained information about the growth rate 
of the instability, and consequently the inhuence of the 


hydrodynamic boundary condition. 


The theoretical approach. The experimental hnd- 
ings can be understood within the lubrication approxi¬ 
mation, from a thin-film model based upon the Laplace 
pressure-driven Stokes equation. We assume incompress¬ 
ible flow of a viscous Newtonian liquid film, of thickness 
varying from 5 to 93 um which is well above the film 
thickness where disjoining pressure plays a role (a few 
tens of nanometers) [34]. Gravitational effects can be 
neglected since all length scales involved in the problem 
are well below the capillary length 4 — 1-73 mm. Finally, 
the velocities of the liquid films are small (e.g. the fastest 
observed rate of change in the amplitude is ~ 25nm/s). 
Thus, the Reynolds and Weissenberg numbers are orders 
of magnitude smaller than 1 and inertial and viscoelastic 
effects can be ignored. 

We non-dimensionalise the problem (see Fig. 1 for vari¬ 
able definitions) through: 
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where the capillary velocity 7/7 is the ratio of the liquid- 
air surface tension to the viscosity of PS. By assuming 
volume conservation, no stress at the liquid-air interface, 
and the Navier slip condition at the solid-liquid bound¬ 
ary, one obtains (see Supplementary Methods) the gov¬ 
erning equation for the dimensionless profile H{X,T): 
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and where the prime denotes the partial derivative with 
respect to X. Note that Eq. (2) is a composite equation 
in the sense that we have kept a second order lubrication 
term in the pressure contribution: the axial curvature. It 
is the lowest order term counterbalancing the driving ra¬ 
dial curvature, and it is thus crucial to obtain the actual 
threshold of the instability. The dynamical aspect of the 
flow is however well described at the lowest lubrication 
order. An interesting discussion on this matter can be 
found in Ref. [9]. 

Performing linear stability analysis, namely letting 
H{X,T) = Hq + e{T)e''^^, where e{T) <C 1, yields 
an exponential growth of the perturbation of the form 
e(T) (X ^ where the rate function is given by: 
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We dehne the fastest growing mode Q* corresponding 
to the smallest time constant t* = t{Q*), and obtain 
Q* = l/(77gv^). Writing Q* in terms of the dimension¬ 
less wavelength leads to: 


A* = 2ttV2Ho , 


( 5 ) 
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FIG. 2. Influence of the geometry on the wavelength of 
the instability. Wavelength of the fastest growing mode as 
a function of the total initial radius (see Fig. 1 (a)). The black 
dashed line represents Eq. (5). The error bars are calculated 
from the error in the geometry and the inaccuracy given by 
the wavelength measurement. 



which is similar to the classical PRI dominant wavelength 
for inviscid jets with no solid core. We note that the 
correspondence between Rayleigh’s result and ours is de¬ 
pendent on both the assumption of Stokes flow and on 
the presence of a solid core in our system. 

The dimensionless growth rate of the fastest growing 
mode is given by: 

= aR -b ^ , (6) 

T* 

where B is the dimensionless slip length, defined in 
Eq. (1), and where a and /3 are parameters that depend 
on geometry only: 

As one sees, the wavelength of the fastest growing mode 
depends exclusively on the initial total radius, while the 
corresponding growth rate is a linear function of the slip 
length. Note that in the case of shear-thickening liq¬ 
uids, the PRI should be significantly slowed down as a 
result of the increasing viscosity associated with an in¬ 
creasing strain rate. Conversely, a shear-thinning liquid 
is expected to accelerate the rise of the instability. In the 
case of a viscoelastic material, Maxwell-like rheological 
models [35] can be implemented and may reveal interest¬ 
ing physics, beyond the scope of the present study. 

The spatial evolution of the instability. Figure 2 
displays the wavelength A* of the fastest growing modes, 
measured on no-slip and slip fibres, as a function of the 
initial total radius hq. As expected from Eq. (5), the 
wavelength A* grows linearly with increasing radius of 
the fibre-polymer system, and is identical on slip and no¬ 
slip fibres. This is consistent with experiments and a the¬ 
oretical framework for retracting liquid ridges on planar 


FIG. 3. Temporal growth of the perturbation. Semi- 
logarithmic plot of the evolution of the perturbation on (a) 
a no-slip fibre (glass), and (b) a slip fibre (glass coated with 
AF2400). The radius of the fibre a and initial thickness of 
the polymer film eo (see Fig. 1 (a)) are indicated. The solid 
line is the best linear fit in the initial regime. 

substrates [26]. The spatial morphology of the instabil¬ 
ity at short times is thus unaffected by the solid-liquid 
boundary condition. While it is clear from Fig. 2 that the 
wavelength is the same on slip and no-slip fibres, there 
is a small systematic deviation from the theory. This 
slight deviation could perhaps be related to the lowest 
lubrication order of the present model [9], but could also 
be attributed to the experimental contribution of sev¬ 
eral modes and the asymmetry of the rate function (see 
Eq. (4)) in the vicinity of the fastest growing mode. 

The temporal evolution of the instability. We 

now turn from the spatial morphology of the instability 
to the temporal evolution. From the experimental im¬ 
ages (see Fig. 1 (b)), we extract the maximal radius of 
an individual bulge as it develops, in order to obtain the 
amplitude C as a function of time. The linear stability 
analysis presented predicts a perturbation that grows ex¬ 
ponentially with a dimensionless growth rate 1/r* for the 
fastest growing mode (see Eq. (6)). Figure 3 displays typ¬ 
ical data for the logarithm of the perturbation amplitude 
normalised by the radius of the fibre. C/a, as a function 
of t, for both no-slip and slip fibres. The data for both 
boundary conditions is consistent with the expected ex¬ 
ponential growth in the early regime. Thus, the initial 
slopes of these curves provide reliable measurements of 
the growth rates. 

The dimensionless growth rates l/r* are shown in 
Fig. 4 for both slip and no-slip fibres, as a function of 
the dimensionless initial total radius Hq. We see that 
for both the slip and no-slip boundary conditions, the 
growth rates show a similar geometry dependence. The 
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FIG. 4. Influence of slip on the growth rate. The inset 
shows the dimensionless growth rate l/r* normalised by the 
no-slip case /?, as a function of a/{j3a), see Eqs. (6) and (7). 
The slip length b is obtained from a best linear fit (dash- 
dotted) to the slip data. The error bars are calculated from 
the error in the geometry and the inaccuracy given by the 
growth rate measurement. The main curve shows the dimen¬ 
sionless growth rate of the fastest growing mode on no-slip 
(glass) and slip (AF2400) fibres, as a function of the dimen¬ 
sionless initial total radius Hq = 1 -(- eo/a (see Fig. 1 (a)). 
Open symbols represent growth rates calculated from Eq. (6) 
using h = 4.0±0.4 fim and the respective experimental geome¬ 
tries. Also shown is the theoretical curve for no-slip (Eq. (6), 
with B — 0). Furthermore, the theoretical curve for slip 
(Eq. (6), with B = 0.3) is plotted as a guide to the eye. 


maxima for the slip and no-slip data can be easily un¬ 
derstood: a decreasing growth rate as Hq converges to 
1 is due to the diminishing thickness of the liquid film, 
and thus to the reduced mobility, while the decreasing 
growth rate for large Hq is due to smaller curvatures and 
thus a smaller driving force of the instability. 

According to Eq. ( 6 ) with B — 0 (no slip), a maxi¬ 
mum in the growth rate is expected for Hq — 1 + cq /a = 
5.15. This prediction is entirely consistent with the data. 
Thus, a coated fibre of a given diameter is maximally 
unstable when the ratio of the film thickness to the fi¬ 
bre radius is about eo/a ~ 4, consistent with an earlier 
theoretical study [5]. The capillary velocity is the only 
adjustable parameter in the no-slip case. We obtained 
7/77 = 294 ± 43 /rm/min, which is in agreement with pre¬ 
vious data [36] adjusted to the temperature used here 
through the William-Landel-Ferry equation [37, 38]. For 
the slip case, there are now two free parameters: the cap¬ 
illary velocity 7/77 and the dimensionless slip length B. 
If we take the value of the capillary velocity to be that 
of the no-slip case, we are left with only one true fitting 
parameter. 

To quantify the slip length in our system the growth 
rates normalised by (3 are plotted as a function of a/{(3a) 
(see inset Fig. (4) and Eqs. ( 6 ) and (7)). As expected. 


the no-slip data is consistent with the theoretical pre¬ 
diction of l/r* = (3 for all geometries. In contrast the 
ratio 1/{t*(3) reveals that the amplihcation due to slip on 
the rise of the instability is more pronounced as a/{(3a) 
increases. For the fibre radii used here, large a/{(3a) 
corresponds to small Hq (see Eq. (7)). The stronger in¬ 
fluence of slip observed for smaller values of Hq is due 
to a non-zero velocity of polymer molecules at the solid- 
liquid interface [31, 33]. For the smallest value of Hq, 
corresponding to eo/a ~ 0.4, the slip-induced amplifica¬ 
tion factor of the growth rate is as large as ~ 4. On the 
contrary, in thick polymer films, the impact of slip is di¬ 
minished. From a best linear ht to the slip-data shown 
in the inset of Fig. (4), the slip length is found to be 
b = 4.0 ± 0.4 fj,m. Obtaining a slip length in the range 
of micrometers is in accordance with former studies on 
the dewetting of entangled polymer films from substrates 
with a fluoropolymer coating [33]. 

Having determined the value of b and knowing the h- 
bre radii, we obtain B for each of the experimental ge¬ 
ometries. Based on the theoretical model (Eq. ( 6 )) corre¬ 
sponding growth rates can be calculated and are shown to 
be in excellent agreement with the experimental data (see 
Fig. (4)). The dimensionless slip length B = b/a ranges 
from 0.25 to 0.37 in our experiments. To guide the eye, a 
typical curve calculated with B = 0.3 is shown in Fig. (4). 
We see that for the fibres with the slip boundary condi¬ 
tion the growth rate is larger than in the no-slip case for 
a given geometry, and the maximum of the growth rate 
1/t* corresponds to a smaller value of Hq. As expected, 
a slippery surface facilitates a higher mobility, and hence 
a faster growth of the instability. This enhanced mass- 
transport also explains the horizontal shift of the maxi¬ 
mum of the growth rate: for a given geometry - and thus 
curvature - the mobility of the slip case is increased in 
comparison to no-slip. Therefore, the maximum of the 
growth rate is shifted to lower values of Hq. 


DISCUSSION 

We report on the Plateau-Rayleigh instability of a vis¬ 
cous liquid polystyrene layer on a solid fibre of radius a, 
when the boundary condition at the solid-liquid interface 
is varied between the classical case of no-slip and the rel¬ 
evant situation of slippage. The wavelength of the fastest 
growing mode on a slip fibre shows a linear dependence 
on the initial total polymer-fibre radius, Hq = a -I- eo, 
and is not affected by the boundary condition, consis¬ 
tent with the lubrication theory developed. For both slip 
and no-slip fibres, we observe an exponential temporal 
growth of the instability at short times, and the respec¬ 
tive growth rates show a qualitatively similar geometry 
dependence. In the case of a slip fibre, the geometry cor¬ 
responding to the maximum of the growth rate 1/t* is 
shifted to a smaller value of ho/a, and the rise of the in¬ 
stability is faster due to the added mobility at the solid- 
liquid interface. The slip induced amplification of the 
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growth rate is significant in a parameter range that is of 
paramount technological relevance. The linear stability 
analysis of the thin film equation developed here is in ex¬ 
cellent agreement with the data, valid for all Newtonian 
liquids, and provides a robust measure of two fundamen¬ 
tal quantities: the capillary velocity 7 / 77 , and the slip 
length b. 


METHODS 
Preparation of fibres 

Glass fibres were prepared by pulling heated glass cap¬ 
illary tubes to final radius in the range 10 < a < 25 /rm 
using a pipette puller (Narishige, PN30). To prepare 
the slip boundary condition, the glass fibres were 
hydrophobised by dip-coating in a 0.5 wt% solution 
of AF2400 (Poly[4,5-difluoro-2,2-bis(trifluoromethyl)- 
1,3-dioxole-co-tetrafluoroethylene]) (Aldrich) in a 
perfluoro-compound solvent (FC72"'"^, Fisher Sci- 
entihc), to form an amorphous fluoropolymer layer 
(Tg ~ 240°C). A dip-coating speed of Imm/s resulted 
in an AF2400 layer with a thickness of 14 ± 1 nm. 
The hydrophobised fibres were annealed in a vacuum 
chamber at 80 °C for 90 minutes to remove excess solvent. 


Preparation of homogeneous PS films 

To prepare homogeneous PS hlms, a concentrated 
solution (35 wt%) of atactic polystyrene (Polymer 
Source Inc.) with a molecular weight of 78 kg/mol and 
low polydispersity = 1.05) was dissolved in 

chloroform (Fisher Scientific). A droplet of the highly 


viscous polymer solution was placed between two glass 
slides, forming a meniscus at the edge of the glass 
slides. We note that the chloroform does not dissolve 
the underlying AF2400 coating on the slip fibres. A 
slip or no-slip fibre could then be placed in the middle 
of the gap between the slides and pulled out of the 
droplet with a constant speed using a motorised linear 
translation stage. By varying the pulling speed vq in the 
range 80 < Uq < 150mm/s, we obtained film thicknesses 
that ranged from cq = 5 to 93 /rm after the solvent had 
evaporated. 


Experimental setup 

The as-prepared samples were placed into a heated 
sample cell to initiate the PRI. With two ~ 0.5 mm 
thick spacers, the coated fibres were suspended above a 
reflective Si wafer (to improve contrast), and placed on 
a microscope hot stage (Linkam). A metal ring in di¬ 
rect contact with the hot stage supported a glass cover 
over the sample and Si wafer (see Supplementary Figure 
1 and Supplementary Methods), thereby ensuring good 
thermal contact and temperature control to within 1°C. 
The surface prohles were analysed from the optical mi¬ 
crographs taken at various times using a custom-made 
edge detection software written in MATLAB. 
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Supplementary Figures 


(a) (b) 



Supplementary Figure 1: Schematic of the sample chamber. Representation of (a) the side view and (b) the top view 
of the experimental setup. On top of a hot stage (A), a Si wafer (B) with two teflon spacers (C) is placed. Freely suspended 
filaments are realized by depositing the fibre (D) across the spacers. The fibre is surrounded by a metal ring (E) with a glass 
cover (F). 
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Supplementary Method 


Fibre thin film equation with slip: 

Here, we provide the main steps of the derivation of the fibre thin film equation (see Eqs. (2) and (3) in the article). 
Let us consider a viscous liquid film coating a cylindrical fibre of constant radius a (see Fig. 1 in the article). The 
axis of the fibre, along the x direction, is chosen as a reference for the radial coordinate r. We assume cylindrical 
invariance of the free-surface profile described by r = h{x,t), at position x and time t. Furthermore, we assume 
an incompressible flow of the viscous Newtonian fluid, where gravity, disjoining pressure and inertia are negligible. 
Defining the local excess pressure p{r, x, t) and the velocity field v(r, x, t), and assuming the lubrication approximation 
to first order, namely that the slopes of the profiles remain small and the velocity is mainly oriented along the fibre 
axis, V = vBx, leads to the two projections of the Stokes equation: 

dxP = V (^drrV + (la) 

drP = 0 , (lb) 

where r] is the dynamical shear viscosity that we assume to be homogeneous and constant. The lubrication approx¬ 
imation is valid since the profile slopes remain small in comparison to 1 at early times - the regime in which the 
linear analysis is performed - and since the typical horizontal length scale of the flow is expected to be larger than the 
minimal wavelength of the Plateau-Rayleigh instability, that is given by Amin = Stt/iq ^ eg. This is indeed always the 
case as even in the worst case geometry, cq « 5a, one has eo/Amin ~ 0.13 ^ 1. According to Supplementary Eq. (lb), 
the excess pressure is invariant in the radial direction, to first order in the lubrication approximation. Thus, the excess 
pressure is set by the Laplace boundary condition at the free surface which, within the small slope approximation, 
reads: 


p = 



( 2 ) 


where 7 denotes the air-liquid surface tension, taken to be homogeneous and constant, and where the prime denotes 
the derivative with respect to x. The two terms on the right hand side of Supplementary Eq. (2) correspond to the 
radial and axial curvatures of the free surface. Note that, as already mentioned in the body of the article, we have 
kept a second order lubrication term in the pressure contribution: the axial curvature —h". This is because it is 
the lowest order term counterbalancing the driving radial curvature 1 /ft., and it is thus crucial to obtain the actual 
threshold of the instability [9]. Regarding the boundary conditions, we assume no shear at the free surface and the 
Navier slip condition at the solid-liquid interface, namely: 





(3a) 

(3b) 


where b denotes the slip length. Integrating Supplementary Eq. (1), together with Supplementary Eqs. (2) and (3), 
yields: 


V 


7 (ft' -k ft^ft'") 
Arjh'^ 



25,2 2 , 

— ft^ + — 2ab 
a 


( 4 ) 


for all r € [a, ft]. Volume conservation requires that: 


dth + 


91 . 

ft 


0, 


( 5 ) 


where we introduced the volume flux per radian: 


Q = 


drrv . 


h 


( 6 ) 
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Combining Supplementary Eq. (5) together with Supplementary Eqs. (4) and (6) yields the fibre thin film equation: 


dfh 


Idrjh [ 


4b 


{h' + I log ( “ ) + (-3 j + 4a^ — 8a& + (4o6 — a^) — 


= 0 . 


(7) 


Finally, introducing the dimensionless variables (see Eq. (1) in the article) into Supplementary Eq. (7) yields the 
dimensionless hbre thin film equation (see Eqs. (2) and (3) in the article). 


Experimental setup and details: 


The as-prepared samples were annealed in ambient atmosphere at 180 °C - well above Tg - which causes the 
polystyrene films to melt and the PRI to develop. The evolution of the PRI was recorded with an inverted optical 
microscope (Olympus BX51). At various times, images were taken using a 5x magnihcation objective and a camera 
with a resolution of 1392 x 1040 pixels (Qlmaging, QlClick), resulting in a pixel size of 1.28 /im. To ensure free 
standing filaments, the coated hbres (D) were placed along two ~ 0.5 mm thick teflon spacers (C) located on a Si 
wafer (B) (see Supplementary Figure 1). The Si wafer was used to improve the optical contrast. A metal ring (E) 
in direct contact with the hot stage supported a glass cover (F) over the sample and Si wafer. The temperature 
within the sample cell is uniform (within 1°C) and convection does not affect the measurements. By using spacers 
(C) made of teflon, we ensured that there is no drainage of the liquid PS film towards the edges (PS does not wet 
teflon). The surface profiles were analysed from the optical micrographs taken at various times, using a custom-made 
edge detection software written in MATLAB. 



